Impurity spin texture at a deconfined quantum critical point 
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The spin texture surrounding a non-magnetic impurity in a quantum antiferromagnet is a sensitive 
probe of the novel physics of a class of quantum phase transitions between a Neel ordered phase and 
a valence bond solid phase in square lattice S — 1/2 antiferromagnets. Using a newly developed 
T = Quantum Monte Carlo technique, we compute this spin texture at these transitions and find 
that it does not obey the universal scaling form expected at a scale invariant quantum critical point. 
We also identify the precise logarithmic form of these scaling violations. Our results are expected 
to yield important clues regarding the probable theory of these unconventional transitions 
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A particularly elegant strategy in the study of strongly 
correlated materials exploits the presence of small con- 
centrations of well-characterized impurities in an other- 
wise pure sample. Each impurity acts more or less in- 
dependently of the others to alter the state of the sys- 
tem around it, and these impurity-induced charge and 
spin textures can then be picked up by nuclear magnetic 
resonance (NMR) or scanning tunneling microscopy ex- 
periments. As these local responses are characteristic 
signatures of the underlying low temperature state, such 
experiments provide a valuable window to the underly- 
ing physics, especially if the state in question has strong 
correlations but no obvious charge or spin order. Q 

Some of these experiments 0] have focused on the ef- 
fects of non-magnetic impurity atoms which give rise to 
a missing-spin defect in strongly correlated Mott insula- 
tors. Due to the uncompensated Berry phase associated 
with such a missing moment [4|, it induces a non- trivial 
pattern of spin density around it, and a direct signature 
of this spin texture can be obtained by analyzing the 
pattern of Knight shifts in NMR experiments. Other ex- 
periments have also studied such effects in cuprate high- 
T c superconductors [H]. These experiments have moti- 
vated several theoretical studies of such physics — these 
include calculations of such impurity effects in antiferro- 
magnets @, superconductors pj, as well as at a quan- 
tum phase transition (QPT) signalling the destruction of 
antiferromagnetism. [8[ 

In this Letter, we use impurities to theoretically probe 
a class of unconventional QPTs between an antiferromag- 
netic phase with long range Neel order and a phase with 
Valence Bond Solid (VBS) order in square lattice 5=1/2 
magnets. Using an extension of the Sandvik-Evertz va- 
lence bond projector loop Quantum Monte Carlo (QMC) 
technique Q that two of us have developed recently [l(| , 
we access the total spin 1/2 doublet ground state of the 
system with a missing-spin defect and compute the spin 
texture induced by this impurity at these Neel- VBS tran- 
sitions. We find that this spin texture does not obey the 
universal scaling form expected to hold at scale-invariant 



quantum critical points. Furthermore, by identifying the 
actual logarithmic form of the impurity scaling viola- 
tions, we argue that the data does not support a first 
order transition. 

Much of our interest in the Neel- VBS transition of 
square lattice antiferromagnets stems from the seminal 
work of Senthil et al. [H who argued that this QPT is 
generically continuous, and admits a natural description 
in terms of 'deconfined' strongly interacting S = 1/2 
spinon excitations rather than the order parameter fields 
of conventional Landau theory — in this sense, it falls out- 
side the well-known Landau classification of phase tran- 
sitions, which predicts a first-order QPT. Numerical ev- 
idence for this theoretical proposal of 'deconfined criti- 
cality' is mixed: while Sandvik, llj, [l2[ and Melko and 
Kaul 13[ see an apparently continuous transition con- 
sistent with deconfined criticality 14| in a microscopic 



S=l/2 spin model, Jiang et al. [15[ provide a detailed 
analysis consistent with the competing scenario of a con- 
ventional weakly-first order transition|16l|. 

Our identification of logarithmic violations of impu- 
rity scaling at these Neel- VBS transitions should be con- 
trasted with the near-perfect scaling collapse we observe 
for the spin texture at a different conventional T = 
critical point between a Neel ordered antiferromagnet 
and a quantum paramagnet without spontaneous VBS or- 
der. These results are expected to yield important clues 
about the correct theory of these unconventional Neel- 
VBS transitions. 

We focus here on two putative realizations [h, 12 1 of 
deconfined criticality corresponding to the Hamiltonians 
T~tjQ2 and Hjqs defined by Sandvik and coworkers: 

n JQ2 = - jJ2 p v -Q J2 p ^ Pkl 

<»i> (ij)(kl) 
n.J Q3 = ~JJ2 P ^ ~Q Yl P *J P kl P rs 
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Here, = 1/4 — Sj • Sj, {ij} refers to a nearest neigh- 
bour (n.n.) bond on the square lattice connecting sites 
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FIG. 1: k dependence of |S z (k)| near the k = and k = Q peaks at the (a) Neel-paramagnet and (b)-(c) Neel-VBS transitions. 
Values of the bulk exponent n were taken from Ref. [8| for the Neel-paramagnet transition and from Ref. Il2j for the Neel-VBS 
transitions. Solid lines in a) are fits to power-law forms obtained by using the value r\ f» 0.44 ±0.02 for the impurity exponent, 
consistent with the estimate in Ref0. 



% and j, and (ij)(kl) ((ij)(kl)(rs)) refer to two (three) 
adjacent parallel n.n. bonds. As a foil of the uncon- 
ventional physics of these JQ models, we also study a 
coupled spin-dimer Hamiltonian Hjj' with antiferromag- 
netic n.n. Heisenberg exchange couplings J for all ver- 
tical bonds, and_ J (J ) for even (odd) columns of hor- 
izontal bonds [17[. These models capture two different 
mechanisms for destabilizing the Neel ordered antiferro- 
magnet: while large values of Q favour a VBS phase in 
the JQ2 and JQ3 models, large values of J drive the 
system to a quantum paramagnetic state that has no 
spontaneous symmetry breaking. 

In order to study the impurity physics at these tran- 
sitions, one needs to access the total spin 1/2 doublet 
ground state of an L x L periodic system with one miss- 
ing site (periodic boundary conditions fix L to be even). 
We have adapted [lfj the valence bond projector loop- 
QMC method to enable an efficient computation of 
the properties of ground states with Stot = 1/2, and S* ot 
fixed from the outset to be either +1/2 or —1/2. This ex- 
tension of the projector algorithm performs as well in the 
Stot — 1/2 sector as the original algorithm does in the sin- 
glet sector. Using this modified algorithm, we study the 
Stot — Stot = 1/2 ground states \G) of periodic systems 
with a missing spin at r — at the Neel-VBS transitions 
in n. JQ 2 (at q c = {Q/J)c/((Q/J)c + 1) » 0.962) and 
7~Ljqz (at q c w 0.603) [l2(], and at the Neel-paramagnet 
transition of U Jf (at (/ /J) c « 1.9096) 

The total S z = 1/2 carried by the ground state spreads 
out throughout the sample to form the impurity-induced 
spin texture <&(r) = (G\S z (r)\G). This texture is ex- 
pected to have a smooth uniform part <& u (r) , and a Neel 
component <&" (r) that alternates in sign between the two 
sublattices of the square lattice. 

If the QPT in question obeys standard scaling theory, 
one expects @, El] q>"(r) = ^f u (r/L) and $™(r) = 
^(i+>7)/2 f n { r /L), where 77 is the bulk anomalous exponent 
associated with the Neel order parameter, and /" and 



/" are the scaling forms for the uniform and alternating 
signals. 

Earlier work [8|] has validated this scaling ansatz for 
the conventional Neel-paramagnet QPT by studying two 
coarse-grained fields (representing the uniform and al- 
ternating signals) obtained from the computed texture 
$(r) by a specific choice of coarse-graining procedure. 
Although straightforward to implement, such a proce- 
dure is somewhat ad-hoc, and depends on the choice of 
coarse-graining prescription. 

Here we finesse this difficulty by noting that the Fourier 
transform S z (k) — J2 r ^K r ) ex P(*k ' r ) (with k = 2nm/L 
and m = (m x ,m y ) with integers m x / y = 0,1, ... L — 
1) is expected to have two peaks, one at k = with 
magnitude constrained to be 1/2, and the second one at 
k = Q = (it, 7r) reflecting the tendency to Neel order. 
The standard scaling ansatz [H, [l8| implies that these 
peaks should satisfy the scaling laws 

Sz(q) = 90 (Lq) for |q| < tt/2 
S z (Q + <i) = L (3 -" )/2 <7q(M) for|q<vr/2 (1) 

The advantage of this new k space formulation is clear: 
one may unambiguously test for scaling by simply exam- 
ining the computed S z (k) for k in the vicinity of k = Q 
and k = 0. In particular, the data for 5 z (q) at the transi- 
tion point computed from samples of varying size L must 
fall on top of each other for |q| tt/2. This is a com- 
pletely unbiased test of scaling as it does not need any 
a priori estimate of the bulk anomalous exponent n for 
the Neel order parameter, nor does it rely on a specific 
coarse-graining procedure to define the scaling compo- 
nents of the texture. 

Our first inkling that standard scaling does not work 
at these Neel-VBS transitions comes from the computed 
values of |5 z (q)| shown in FigQ]and Fig[2]for q = 2nm/L 
with |m| <C L/2 (in practice, we focus on |m| < L m ; n /12 
and average over all m that correspond to a given |m|). 
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FIG. 2: 1/L dependence of |5" Z (Q + 2nm/L)\/L (3 ~ v)/2 and \S z {2wm/L)\ for small |m| at the Neel-VBS and Neel-Paramagnet 
transitions. All horizontal lines are guides to the eye that indicate the expected behaviour if scaling was perfectly obeyed, and 
values of r\ are the same as in Fig[T] 



Larger values of L are seen to yield a systematically larger 
value of \S Z \ at the same |m|. This behaviour at the Neel- 
VBS transitions is in clear violation of the scaling form 
Eq. [U this should be contrasted with the excellent scal- 
ing observed at the conventional Nccl-paramagnet criti- 
cal point of the J — J model. Given the unbiased nature 
of this test of scaling, we consider this rather strong evi- 
dence for violation of impurity scaling properties at these 
Neel-VBS transitions in the JQ2 and JQ3 models. 

Next, we analyze the Bragg peak at the antiferro- 
magnetic wavevector, k = Q, focusing on the L depen- 
dence at the Neel-VBS transition point in both JQ mod- 
els. We find no evidence of any double-peak structure 
for the histogram of |S 2 (Q)| at these Neel-VBS tran- 
sitions, implying that the first order jump in the or- 
der parameter, if any, is immeasurably small even at 
sizes as large as L — 96. Furthermore, we confirm 
that the computed values obey the power-law scaling 
\S Z (Q)\ ~ L^-v)/ 2 quite well at both Neel-VBS tran- 
sitions, with the anomalous exponents rjjQ 3 w 0.33 and 



VJQ2 ~ 0-35 taken from Ref 13, 13. Our results for the 
J — J model are also consistent with the power-law scal- 
ing \S Z (Q)\ ~ £( 3 -'/)/ 2 with the known value of 77 « 0.04 
for the Neel-paramagnet QPT Q. 

However, violations of impurity scaling in the stag- 
gered component of the texture at the Neel-VBS tran- 
sitions become evident when one tests for scaling col- 
lapse at k = Q + 2-nva.jL with |m| small but non-zero 



(we focus on |m| < I/ m j n /12 and average over all m that 
correspond to a given |m|) — again, larger L give larger 
values of |5 z (k)| for the same non-zero |m|(Fig Q] and 
Fig©- This is again underlined by the excellent scaling 
collapse exhibited by the corresponding quantities com- 
puted at the Neel-paramagnet quantum critical point of 
the J — J model (Figures Q] and [J) . 
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FIG. 3: Logarithmically modified scaling collapse of S z (k) 
near k = at the Neel-VBS transition, with Iq = 5 ± 1 (lo = 
12 ± 1) for the JQ 2 (JQ3) model. 

What do these violations of quantum-critical scaling 
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FIG. 4: Logarithmically modified scaling collapse of S z (k) 
near k = Q at the Neel-VBS transition, with Iq = 0.75 ± 0.2 
(Iq = 1.5 ± 0.5) for the JQ 2 (JQs) model. 



transition is continuous, the theory of deconfined critical- 
ity [l[ needs to be modified in order to account for these 
logarithmic corrections. This conclusion underscores the 
utility of impurity physics as a probe of complex strongly- 
correlated states of many-body systems. An interesting 
follow-up would be to use the same probe at non-zero 
temperature above the QPT and test for violations of 
scaling predictions for the impurity susceptibility. 
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for a critical reading of this draft. We acknowledge com- 
putational resources of TIFR and funding from DST- 
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signal? This is best addressed by asking if the computed 
spin texture satisfies some suitably modified scaling laws. 
To explore this, we first note that the absence of any 
power-law prefactor to go in the scaling ansatz Eq. [1] re- 
flects the conservation of total spin and the specific power 
of L that multiplies <?q reflects the presence of power-law 
Neel order at criticality, while the scaling argument Lq of 
both functions follows simply from the statement that L 
is the only length scale of relevance to the long-distance 
physics at a scale invariant critical point. 

This form of the scaling argument can break down 
if the effective low-energy theory has a term which is 
marginally irrelevant for the long-distance physics, and 
rcnormalizes to zero slowly. This can introduce an addi- 
tional length scale, and give rise to logarithmic violations 
of scaling. Indeed, signatures of such logarithmic drifts 
have been recently seen by Sandvik in his analysis of the 
bulk Neel-VBS transitions 0. 

This motivates us to ask if the computed spin texture 
obeys a modified scaling form (for |q| <C 7t/2) 

S z (q) = g Q (Lq/\og(L/l a )) 
5 z (Q + q) = L^/ 2 g Q (Lq/log(L/l Q )) (2) 

where lo and Iq are related to the additional non- 
universal length scale introduced by the slow vanishing 
of a marginally irrelevant term in the effective Hamilto- 
nian. As is clear from Figs the answer is yes: this 
modified scaling law gives an extremely good account of 
our results. 

These logarithmic violations of scaling are at odds with 
predictions of the theory of deconfined criticality [lij . 
However, given the absence of any clear signal for phase- 
coexistence at the transition point, one cannot simply 
ascribe these violations to the presence of an immea- 
surably weak first-order jump in the order parameters. 
We are thus led to conclude that although the Neel-VBS 
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